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I. INTRODUCTION 



The XY model in 3 dimensions is an interesting system in statistical mechanics. Physi- 
cally it describes the critical behaviour of superfTuid He4 From the theoretical point of 
view it provides a relatively simple example of an order- disorder phase transition in which 
the condensation of solitons plays an essential role [2-5]. It is fairly well established that the 
transition is continuous, and the basic critical indices are known with good precision . 

A Ginzburg-Landau kind of description of the system has been given, describing phe- 
nomenologically the condensation of vortices in the high temperature phase ||. 

In this paper we shall study the system from a slightly different point of view. We 
shall look at it as the lattice formulation of a 2 + 1 dimensional quantum field theory. We 
shall then show that the (2 dimensional) ground state of the theory spontaneously breaks the 
conservation of the number of vortices at high temperature; vortices condense in the vacuum 
in the same way as Cooper pairs do in a superconductor. In that phase the number of vortices 
is not defined in the vacuum, which is a superposition of states with different vorticity. We 
shall directly exhibit this condensation by showing that the vacuum expectation value (vev) 
of an operator which changes the number of vortices is non zero. 

An analogous technique has been used to detect monopole condensation in the ground 
state of U(l) gauge theory [[Tof and in SU(2) gauge theory, as a mechanism for confinement 
of colour . 

In sect. 2 we shall briefly introduce the model and fix the notation. We shall also define 
a creation operator for a vortex. 

In sect. 3 we shall present and discuss the results of numerical simulations where the vev 
of the above operator is measured. The result will be a direct detection of the spontaneous 
breaking of the conservation of the number of vortices in the disordered phase and a direct 
evidence of vortex condensation. We shall also obtain an alternative determination of known 
critical indices, which agrees with existing results, and measure the critical index related to 
our disorder parameter. 



II. THE CREATION OPERATOR OF VORTICES. 

The lattice action of the system is 

S = /?]T(l-cos(A^))) (1) 

fj,,i 

The field variable is the angle 9 in the site i: /i runs from to 2, being the time axis. 
In the limit of zero lattice spacing (a — > 0) 

s-f E(V0V + cV) (2) 

ft,i 

At high /3's the system describes a massless free scalar particle. The couplings of higher 
dimension appearing in the series expansion of the cosine become important at lower (3. 
At (3 C ~ 0.454 the system undergoes a phase transition, which is known to be of second 
order ||. 
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The system has soliton configurations with the geometry of vortices: numerical simula- 
tions show that vortices play an important role in the phase transition 
A vector field 

= (3) 

can be defined, and a current = e m pd a Ap associated with it, which is trivially conserved 
in smooth configurations: 

d^j, = Ve^dTA? = (4) 

because of the antisymmetry of the Ricci tensor. Eq(|]) is the analog of the Bianchi identity 
in QED. 

The constant of the motion is 

$ = J d 2 x j°(x, x°) = J d 2 x (V A A) (5) 
On the other hand from eq. (|3|) 

$ = j> A-dx = n-2n (6) 

with C any closed path. 

A vortex is defined as a configuration for which n is non zero. 

0(x - y) = arctg^ — ^ (7) 

{x-y)i 

is an example of such a configuration with n = 1. It has a singularity at x = y. 

Eq.'s tell that the number of vortices is a constant of motion, and define a U{1) 
symmetry. 

We shall show that, for (3 > (3 C this symmetry is realized a la Wigner: for (3 < (3 C a 
spontaneous breaking of it occurs. 



To do that we will define an operator which creates a vortex In the continuum 

version of the theory 8q9 is the conjugate momentum to the field variable 9. To add a vortex 
in the site y, 6(x — y) to any field configuration the analog of the simple translation 



x) = \x + a) 



can be used. The translation operator 

n(x, t) = exp 



(8) 



6 yd 6(y,t)6(x-y) 

indeed, when applied to the state \9(x, t)) gives 

fx(x,t)\9(y,t)) = \9(y,t)+9(x-y)) (9) 
The operator number of vortices 
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V(t)= I d 2 xf(x,t) 



(10) 



has commutation rule with \i 

\V(t),fi(x,t)] = I d 2 z V A A(z, i) /i(x, t) = 2vr/i(f,t) 



(11) 



where A = V0. 

If the U(l) symmetry which counts the vortices is realized a la Wigner, then the ground 
state must have a definite number of vortices and (/x) = 0. Since /i changes the number of 
vortices a vev (//) ^ signals their condensation. 

The naive translation of eq. (§) on the lattice is [T(| 



n(n, no) = exp 



-(3 sin(A o 0(n*, n Q ))9(n" - n, n Q ) 



(12) 



where the site n' runs on the sites of a 2 dimensional slice of the lattice at constant no, 
except the position of the vortex. 



Instead of Eq.(12) we shall use a compactified version of it. By this we mean 



pi = exp 



(cos(A 6 + 6)- cos(A 6 



(13) 



which coincides with eq.(^) at first order in 9. With the definition ( pT3[) (/i) is independent 
of the choice of the gauge for A^, i.e. of the zero for the angle 9: a redefinition of it is 
reabsorbed in the integration variable. 
Two remarks are in order 

1) in principle one expects (/i) = in the ordered phase (3 > (3 C . In fact (/i) is an analytic 
function of j3 if the number of degrees of freedom is finite, and therefore, if it were zero 
for /3 > (3 C it would be identically zero for all /3's. Only in the limit of infinite volume 
(thermodynamic limit) singularities develop, and (fi) can be identically zero above /3 C , 
without being zero everywhere. Indeed a simple computation at large f3's, where the 
integral is approximately gaussian, gives 



(ji) ~ exp 



-P (ciV^ + ^ + O^ 



(14) 



2) Due to periodic boundary conditions the total number of vortices must be zero as 
easily seen from eq.@ with the boundary C at spatial infinity. In principle then the 
correct way to determine (/i) is to observe a vortex antivortex pair correlation 

(M^MoXcJ^ 2 ( 15 ) 

By cluster property and C invariance it tends exponentially to (fi) 2 at large distances. 
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In practice, since in the thermodynamical limit physics becomes independent of the 
boundary conditions and decorrelated from them at distances larger than the corre- 
lation length, what happens is that also a single vortex can be created, and a direct 
measurement of (p) can be done. (p) determined in this way is indeed independent 
of the boundary conditions within errors. When a single vortex is put on the lattice 
by our operator a dislocation appears at the boundary with opposite vorticity to cor- 
rect for the boundary conditions, and, if the lattice is big compared to the correlation 
length, what is obtained is the same (/i) 2 as with a vortex antivortex pair at large 
relative distance. 



III. NUMERICAL SIMULATION AND RESULTS. 



Simulations have been done on a 20 2 x 40 and 30 2 x 60 lattices to measure the (pp) 
correlation. Single vortices (p) have also been measured on the same lattices and on 20 3 , 30 3 
and 40 3 lattices. Instead of directly determining (p) as a function of (3 it is more convenient 
to measure JTT) 



and to reconstruct (p) 2 as 



(16) 



cxp 



pifi) W 



(17) 



The observed shape of p is shown in fig.l. The sharp negative peak at (3 C indicates a steep 
drop of (p) toward zero at (3 C . At low /3's p is compatible with and independent of V 
(fig.2). This means by eq.([H]) that limy_ >00 (/i) ~ 1 for (3 < (3 C , and thus that vortices do 
condense in the vacuum. At (3 > (3 C the behaviour of eq.(|l4D can be checked by numerical 
calculation of the gaussian integral. For a cubic lattice with a monopole in the center the 
result is 



p= -11.332 -L + 72.669 + 0{l//3) (18) 

The comparison with numerical simulations is shown in fig.l. 

Eq . fll~8|) means that (p) 2 is exactly zero in the limit V — > oo. Eq.(|T8D agrees with the 
observed values of p at large /3's (fig.l). 

Near f3 c where the correlation length goes large, a finite size analysis can be done to 
explore the limit of infinite volume and determine the critical index 5 of p, /3 C , and the 
critical index v of the correlation length. 

As ((3 - (3 C ) -»• 0~ and V -> oo 

(p) ~ ((3 C - P) 5 (19) 

At finite L, 
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^w|;j (20) 

As (3 — > /3 C £ diverges, £ ~ (/3 C — /3) _l/ and the dependence on a/£ can be neglected 

(ri^f(L^(p c -P)) (21) 



or 



p^L^f^'iPc-P)) (22) 

p/L l l u in the scaling region must be a universal function of L 1//!/ (/3 c — /3). The quality of this 
scaling is shown in fig. 3. The critical index v and the critical temperature (3 C are determined 
by enforcing this scaling behaviour. 

Fig. 4 shows the data for a vortex antivortex pair at large distance, compared to the 
single vortex data. In the condensed phase the single vortex behaves in fact as a pair at 
large distance, confirming that there is no problem in the thermodynamical limit. At (3 > (3 C 
the two quantities are different and at f3 sufficiently large are well described by eq. flll]) with 
the different appropriate values of C\ and C2, as shown in the figure. 

A more precise finite size scaling analysis assuming a specific form for / for (3 < (3 C 

f{L l/v {(3 c ~I3))-- " 



(L l /»((3 c -f3)+v 1 ) 2 + v 2 



gives a precise determination of the critical indices and of f3 c 

5 = 0.740 ± 0.029 (24) 
(3 C = 0.4538 ± 0.0003 (25) 
v = 0.669 ± 0.065 (26) 

with x 2 /d°f — 1-07. For comparison the determination of ref. [8] is v = .670(7) and (3 C = 
.45419(2). 

The behaviour of the correlator (p(x)fi(O)) as a function of the distance gives an estimate 
of the correlation length ^m{(3) 

(ji(x)n{0)) - (fi) 2 ~ exp [-x/£ M (J3)] (27) 

£m is of the order of a few lattice spacings (see fig.5). £m(P) should scale as {(3 C — [3)~ u as 
((3 — (3 C ) — > _ . With our relatively large statistical errors we were not able to check this 
behaviour. 



IV. CONCLUDING REMARKS. 

We have constructed a disorder parameter for the A transition in the 3c? XY model which 
is the vev of an operator with non trivial vortex number. It is different from zero at high 
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temperatures (ft < ft c ) in the disordered phase, thus demonstrating condensation of vortices. 
In the ordered phase it tends to zero exponentially with the linear size of the lattice. 

ft c and the critical index v of the correlation length can be determined with good precision 
by a finite size scaling analysis of our disorder parameter (//). As (ft — ft c ) — * 0~ /x ~ (ft c — ft) 5 
with 5 = .74 ± .03. 

A similar construction also works for other models undergoing an order disorder transi- 
tion produced by condensation of solitons. An example is the condensation of monopoles in 
gauge theories [fT0| , |Tlj| . 
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Figure Captions 

Fig.l p as a function of (3. The dashed lines at high values of j3 are the perturbative estimates. 

Fig. 2 p as a function of 1/L (L = size of the lattice) in the condensed phase, p converges 
to a value consistent with zero. The abscissae for different betas are slightly displaced 
for graphycal reasons. 

Fig. 3 Quality of the finite size scaling analysis. The figure corresponds to optimal values of 
j3 c and v. 

Fig. 4 Comparison of p for single vortex to p for vortex - antivortex pair at large distance. 
The two curves coincide within errors in the condensed phase. They differ at high f3 
and agree with the perturbative determination (dashed lines on the right side). 

Fig. 5 p as a function of distance at two values of (3. The estimated correlation length is few 
lattice spacings. 
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